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Abstract 

The question as to what the relevant effective degrees of freedom at the chiral phase tran- 
sition are remains largely unanswered and must be addressed in confronting both terrestrial 
and space laboratory observations purporting to probe matter under extreme conditions. 
We address this question in terms of the vector susceptibility \v (VSUS in short) and the 
axial-vector susceptibility xa (ASUS in short) at the temperature-induced chiral transition. 
We consider two possible, albeit simplified, cases that are contrasting, one that is given by 
the standard chiral theory where only the pions figure in the vicinity of the transition and the 
other that is described by hidden local symmetry (HLS) theory with the Harada-Yamawaki 
vector manifestation (VM) where nearly massless vector mesons also enter. We find that 
while in the standard chiral theory, the pion velocity proportional to the ratio of the 
space component of the pion decay constant over the time component /* tends to zero 
near chiral restoration with /* ^ 0, in the presence of the vector mesons with vanishing mass, 
the result is drastically different: HLS with VM predicts that xv automatically equals xa 
in consistency with chiral invariance and that u„ ~ 1 with /* w /* — > as T — > T c . These 
results are obtained in the leading order in power counting but we expect their qualitative 
features to remain valid more generally in the chiral limit thanks to the VM point. 

*Present address: Department of Physics, Nagoya University, Nagoya 464-8602, Japan. 

'and Department of Physics and Astronomy, University of South Carolina, Columbia, SC 29208 



1 Introduction 



One of the most crucial questions to answer in the effort to understand chiral restoration in 
relativistic heavy-ion collisions as well as in dense medium as in compact stars is: What are the 
relevant effective degrees of freedom just before and after the phase transition? The standard 
scenario, generally accepted in the community, is that the only relevant excitations in the broken 
symmetry sector near the phase transition are the pions, the pseudo-Nambu-Goldstone modes 
of broken chiral symmetry, i.e., the standard chiral theory. However there is no a priori reason 
to exclude other scenarios. In fact, it has been argued by Harada and Yamawaki p^E] that the 
vector manifestation (VM) with the massless (in the chiral limit) vector mesons can correctly 
describe chiral restoration. In a recent attempt to understand some of the puzzling results 
coming out of relativistic heavy ion experiments at RHIC, Brown and Rho [3] invoked the VM 
scenario in which the "light" p mesons play a crucial role: The vector mesons there are "relayed" 
via a Higgsing to the gluons in the QCD sector. 

In this paper, we address the issue of what the relevant degrees of freedom can be at the 
chiral transition induced by high temperature and their possible implications on observables in 
heavy- ion physics. In doing this, we focus on the vector and axial- vector susceptibilities and 
the pion velocity very near the critical temperature T c using the result obtained in |4j who have 
shown that the VM holds at T = T c . The issue of what happens at high density is discussed in 



As a way of introduction to the main objective of this paper, we begin by briefly summarizing 
the arguments and results obtained in the standard chiral theory scenario Ej ■ 

Consider the vector isospin susceptibility (VSUS) xv (denoted by SS as \i) an d the axial- 
vector isospin susceptibility (ASUS) \A (denoted by SS as X15) defined in terms of the vector 
charge density Vq(x) and the axial- vector charge density Aq{x) by the Euclidean correlators: 




(1.1) 



(1.2) 



where ( )p denotes thermal average and 



V a = ^° Y ^, A\ 




(1.3) 



with the quark field ip and the r a Pauli matrix the generator of the flavor SU(2). 
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We are interested in these SUS's near the critical temperature T = T c at zero baryon density 
n = 0. In particular we would like to compute them "bottom-up" approaching T c from below. 
In order to do this, we need to resort to effective field theory of QCD which requires identifying, 
in the premise of an EFT, all the relevant degrees of freedom. 

Let us first assume as done by Son and Stephanov (SS) [HI E] that the only relevant effective 

degrees of freedom in heat bath are the pions, and that all other degrees of freedom can be 

integrated out with their effects incorporated into the coefficients of higher order terms in the 

effective Lagrangian. Here the basic assumption is that near chiral restoration, there is no 

instability in the channel of the degrees of freedom that have been integrated out. In this pion- 

only case, the appropriate effective Lagrangian for the axial correlators is the in-medium chiral 

Lagrangian dominated by the current algebra terms, 

ft ^ i 
Leff = *f- (TWot/Vo^ - vlTrdiUdirf) - -{^)ReM^U + ■■■ (1.4) 

where is the pion velocity, M is the mass matrix introduced as an external field, U is the 
chiral field and the covariant derivative VoU is given by VqU = OqU — ^ha(t~3U + Ut%) with 
\jla the axial isospin chemical potential. The ellipsis stands for higher order terms in spatial 
derivatives and covariant derivatives. ^ Given the effective action described by (|1.4j) with 
possible non-local terms ignored, then the ASUS takes the simple form 

d 2 . t 2 -s 

XA - —Q^Wff\^ A =o = U ■ (1-5) 

The principal point to note here is that as long as the effective action is given by local terms (sub- 
sumed in the ellipsis) involving the U field, this is the whole story: There is no other contribution 
to the ASUS than the temporal component of the pion decay constant. 

Next one assumes that at the chiral phase transition point T = T c , the restoration of chiral 
symmetry dictates the equality 

XA = XV- (1-6) 

While there is no lattice information on xa, XV has been measured as a function of tempera- 
ture [HUH]. In particular, it is established that 

* lr The notation here deviates a bit from that of SS. For example, it will turn out that the pion velocity will 
have the form v% — f%/ft ( see Eq. 15.14H ') where /* (/£) is the temporal (spatial) component of the pion decay 
constant. 
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Xv\t=t c + 0, (1.7) 
which leads to the conclusion [7j that 

fl\T=T c + 0. (1.8) 

On the other hand, it is expected and verified by lattice simulations that the space component 
of the pion decay constant should vanish at T = T c . One therefore arrives at 

«2~/*/#-0, T^T C . (1.9) 

This is the main conclusion of the standard chiral theory. 

To check whether this prediction is firm, let us see what one obtains for the VSUS in the same 
effective field theory approach. The effective Lagrangian for calculating the vector correlators 
is of the same form as the ASUS, Eq. (|1.4j) . except that the covariant derivative is now defined 
with the vector isospin chemical potential [iy as VqU = OqU — \^v{ T ?,U — Ut^). Now if one 
assumes as done above for \A that possible non-local terms can be dropped, then the SUS is 
given by 

d 2 

XV = ~gJ^T Le ff\^v=0 I 1 - 10 ) 
which can be easily evaluated from the Lagrangian. One finds that 

XV =0 (1.11) 

for all temperature. While it is expected to be zero at T = 0, the vanishing \v f° r T ^ is at 
variance with the lattice data at T = T c . # 2 

We now turn to the main objective of this paper: the prediction by the vector manifestation 
(VM) 1 . Basically the same scenario was suggested some time ago in conjunction with Brown- 
Rho scaling jlUl lllj . As will be shown in detail in the following sections, the VM requires that 
the vector mesons figure on the same footing with the pions as the relevant degrees of freedom as 
the chiral transition point is approached from below. The key reason for this conclusion is that 
the chiral transition coincides with the VM fixed point at which the vector meson mass must 
vanish in the chiral limit |12j . This means that the vector-meson degrees of freedom cannot be 
integrated out near chiral restoration. 



#2 



The reason for this defect is explained in terms of hydrodynamics by Son and Stephanov 
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Our principal results - which are basically different from the standard chiral theory scenario 
- can be summarized as follows. In the presence of the p-meson, the only approach that is 
consistent with chiral perturbation theory is the hidden local symmetry (HLS) with the VM 
fixed point * 3 . The present analysis based on this theory predicts # 4 

fi\T=T c = ti\T=T c = , V n \ T=Tc < 1 (1.12) 

and 

N f 2 

Xa\t=t c = Xv\t=t c = -£-T c , (1.13) 

O 

where we have included the normalization factor of 2Nf. Note that the equality of xa and xv 
at T = T c is an output of the theory. This result is a direct consequence of the fact that the p 
and 7r enter on the same footing in the VM: At the VM fixed point, the longitudinal components 
of the vector mesons and the pions form a degenerate multiplet. 

The rest of the paper is devoted to the derivation of the main results (fl~T2l) and fTHty . In 
Section hidden local symmetry (HLS) theory is briefly introduced. Section |3] describes how 
thermal two-point functions are calculated in the HLS theory. In Section^ we write down the 
in-medium vector and axial-vector current correlators that are needed in what follows. Pion 
decay constants and pion velocity are computed in the given framework in Section |SJ The 
susceptibilities are defined in Section and computed for temperature T ~ T c . The conclusion 
is given in Section The Appendices contain explicit formulas used in the main text. A more 
extensive treatment of the material covered in this paper together with other issues of finite 
temperature effective field theory in the VM is found in Ref. |14j . 



2 Hidden Local Symmetry 

In this Section, we briefly summarize the HLS model. Our discussion will be highly sketchy. 
For details, the readers are invited to the review j^j- As mentioned above, the HLS spin-1 field 

is assumed to be as relevant as the pion field near chiral restoration. The HLS model |15| I16j 

* 3 It has been stressed in the literature (see, e.g., 1 131 [5]) _ and is stressed again - that HLS is a bona-fide 

effective field theory of QCD only if the p-meson mass is considered as of the same chiral order as the pion mass. 

In HLS theory, this condition is naturally met by the p-meson mass near the chiral transition point, so chiral 

perturbation theory should be more effective in this regime. This point that underlines our arguments that follow 

justifies our one-loop calculation. 
# 4 The reason for that the v n deviates from 1 will be explained below. 
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is based on the G g i bai x -f^iocai symmetry, where G = SU(iV/) L x SU(7V/) R is the global chiral 
symmetry and H = SU(iV/-) v is the HLS. The basic quantities are the HLS gauge field and 
two variables or "coordinates" 

a,R-^e^, (2.1) 

where tt denotes the pseudoscalar Nambu-Goldstone (NG) boson and a the NG boson absorbed 
into the HLS gauge field (longitudinal p). F n and F a are corresponding decay constants, and 
the parameter a is defined as a = F%/F%. The transformation property of £l,r is given by 

&,r(s) - CL, R (x) = fe(x)e L ,R(x)ri )R , (2.2) 

where /i(x) 6 -Hiocai and <7l,r £ ^global- The covariant derivatives of £l,r are defined by 

D^L = - iV^L + i£,h£>n , 

D^r = d^ R - iV^ R + , (2.3) 

where £ M and 7£ M denote the external gauge fields gauging the Global symmetry. From the 
above covariant derivatives two 1-forms are constructed as 

«l = (^R-e R -^L-ei)/(^) , 

«f = (D^ R -^ + D^ L -^)/(2i) . (2.4) 



It should be noticed that, as first pointed by Georgi in Ref. [12] and developed further in 
Refs. [171 1181 12*]. the systematic chiral perturbation can be performed with including the vector 
meson loop in addition to the pion loop in the HLS. The expansion parameter is a ratio of the 
p meson mass to the chiral symmetry breaking scale A x , m p /A x , in addition to the ratio of the 
momentum p to A x , p/A x , as used in the ordinary chiral perturbation theory. The counting 
scheme is made as in the ordinary chiral perturbation theory by assigning 0(p) to the HLS 
gauge coupling g ^31 El ■ 

With the above counting scheme the Lagrangian at the leading order, counted as 0(p 2 ), is 
given by ^3 

C = Fl tr [a^] + F 2 a tr [d^df 1 - -Ltr [V^V^] , (2.5) 

Zg 

where g is the HLS gauge coupling and = d^V u — d u V^ — i[V^, V u ] the gauge field strength. 
When the kinetic term of the gauge field is ignored in the low-energy region, the second term of 



Ea. H2.5JI vanishes by integrating out and only the first term remains. Then, the HLS model 
is reduced to the nonlinear sigma model based on G/H. 

The one-loop quantum corrections calculated from the leading order Lagrangian in Eq. IJ2.5JI 
are counted as 0(p ). The divergences generated at 0(p 4 ) are renormalized by the 0(p 4 ) terms, 
a complete list of which are given in Refs. |17l l2"]. Here we show the terms of the 0{p A ) Lagrangian 
relevant to the present analysis J171 118j : 

£ (4) = zM%vVn + Z2tx[A tiv A^} + z 3 ti[V^Vn, (2-6) 

where 

v»v = \ [i R n^ R + SlC^sI] , (2.7) 

with 1Z^ U and C^ u being the field strengths of IZ^ and C^. 

We should stress that we assume that we obtain the bare HLS Lagrangian by integrating out 
the quark and gluon degrees of freedom at the matching scale A, so that bare parameters of the 
above HLS Lagrangian such as F w bare defined at A are determined through the Wilsonian match- 
ing between the HLS and the underlying QCD ^H]: As we briefly review in Appendix lA.il bare 
parameters are determined through the Wilsonian matching conditions (JA.6JI - IJA.8JI obtained by 
matching the axial-vector and vector current correlators in the HLS with those in the operator 
product expansion (OPE). As was shown in Refs. |18l l2j. these bare parameters are scaled down 
to the low energy region through the Wilsonian renormalization group equations (RGEs) to 
predict several physical quantities in remarkable agreement with experiment. 

Let us extend the above HLS Lagrangian to the analysis in hot matter. We assume that we 
obtain the bare HLS Lagrangian by integrating out the quarks and gluons degrees of freedom at 
the matching scale A in the presence of medium, and then the bare parameters are determined 
by matching the HLS to the underlying QCD. As was shown in Ref. |I] and briefly reviewed in 
Appendix IA.21 when we make the matching in the presence of hot matter, the bare parameters 
have the intrinsic temperature dependences. In general, the Lorentz non-scalar operators such 
as qj^D^q exist in the form of the current correlators derived by the OPE However, 
as discussed in Appendix IA.2I such Lorentz symmetry violating contributions caused from the 
Lorentz non-scalar operators are suppressed by, at least, a factor of 1/A 6 compared with l+a s /ir, 
and the Lorentz violating effects in the bare ir decay constant and the bare a (longitudinal p) 
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decay constant are small: At bare level, the difference between i 7 * bare and bare as well as that 
between F* h&re and F^ haxe is small. Furthermore, since we will study the physical pion decay 
constants and the vector and axial-vector susceptibilities only near the critical temperature, 
the Lorentz violating effect in the HLS gauge coupling, which may distinguish from gi (see 
Appendix IA.2j) . is irrelevant due to the decoupling nature of the transverse p near the critical 
temperature in the VM. Thus, here we use the Lagrangian Q2.5JI with Lorentz invariance even at 
non-zero temperature to calculate the quantum and hadronic thermal corrections. We should 
stress that the explicit Lorentz violation in medium - which is not negligible - is of course taken 
into account (see Appendices). 



3 Two-Point Functions 

We need to consider two-point functions involving the isovector vector and axial-vector cur- 
rents. Note that, when we calculate the hadronic thermal corrections, we assigned 0{p) to the 
temperature T as in the approach based on the ordinary chiral perturbation theory (20]. * 5 
Let us calculate them to one- loop order as in Refs. ^3 |2]. In calculating the loops, we adopt 
the background field gauge (see Refs. |18| 12] for details in HLS theory) and the imaginary time 
formalism (see, e.g., Ref. (23)- For convenience, we introduce the following Feynman integrals 
to calculate the one-loop hadronic thermal corrections and quantum corrections to the two-point 
functions: 

°° r d 3 k 1 
d 3 k 



n=—oo 
oo 



f d k 1 

B ip ,p;M u M 2 ;T)^Tj:j (27r)3 [M * - k*\[M* - (k - p)*] ' ^ 



n=— oo 
oo 



u^i - m m t\-t sr ( d3k {lk-pf{2k-py 

(Po,P;MuM 2 ;T)=Tj:j (27r)3 [M? _ fe2] [M | _ {k _ p?] > ( 3 - 3 ) 

where p = \p\ , and the Oth component of the loop momentum is taken as A; = i2mrT, while that 
of the external momentum is taken as p° = i2n'irT [n'\ integer). Using the standard formula 
(see, e.g., Ref. [221) > we can convert the Matsubara frequency sum into an integral over ko with 
ko taken as the zeroth component of a Minkowski four vector. Accordingly, the above functions 
are divided into two parts as 

A (M; T) = ^(M) + A (M; T) , 



#5 The same treatment within the framework of the HLS was done before in Refs. |21 1 P|. 




Figure 1: Feynman diagrams contributing to the A^-Ay two-point function. Here n represents 
the quantum pion field and likewise for the others. 

B (p ,p;M 1 ,M 2 ;T) = B^ (p ,p; M t , M 2 ) + B (p ,p; M u M 2 ;T) , 
B^(p ,p;M 1 ,M 2 ;T) = B^"" (p ,p; M Xl M 2 ) + B^(p ,p; M u M 2 ;T) , (3.4) 

where ^ vac) , B^ c) and #( vac )^ are given by replacing T E^-oo with J §t in Ec l s - CEIMII3J> 
and Aq, Bq and B^ u are defined by Eq. (|3.4jl . In the present analysis, the forms of i4.jj , 
Bq and B^ vac> ^ 1 ' are equivalent to the zero-temperature ones. Then, with po taken as the 
Oth component of the Minkowski four vector, they have no explicit temperature dependence 
while the intrinsic dependence mentioned above remains. Therefore, the functions ^4Q Vac \ Bq 7 ^ 
and j g( vac )M 1/ represent quantum corrections. In Bq and B^ v one can perform the analytic 
continuation of po to the Minkowski variable after integrating over k^: Here po is understood as 
Po + ie (e — > +0) for the retarded function and po — ie for the advanced function. It should be 
noticed that the argument T in the above functions refers to only the temperature dependence 
arising from the hadronic thermal effects and not to the intrinsic thermal effects included in the 
parameters of the Lagrangian. 

Now, let us calculate the one-loop corrections to the two-point function of the axial-vector 
(background) field A^. This is obtained by the sum of one particle irreducible diagrams with 
two legs of the axial- vector background field A^ . In Fig. ^ we show the Feynman diagrams 
contributing to the A^-Ay two-point function at one-loop level. With the help of 1)3. 4 J) , one can 
express the one-loop corrections to the two-point function in a simple form. The corrections 
from p and/or tt shown in Figs.^a)-(c) lead to the two-point function 

U { l- loop) ^(p ,p;T) = -N f aM 2 p g^B (p ,p;M p ,0;T) + N f ~B^(p ,p; M p ,0;T) 

+ N f (a-l)g^A (0,T) , (3.5) 
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where i?o-term comes from Fig. ^a), and S^-term and Ao-term from Fig. ^b) and Fig. ^c), 
respectively. Combining the above loop corrections with the tree contribution given by 

nf ee) ^( P0 ,p) = <T F 2 bare + 2z 2 , barc (^V - V »f) , (3.6) 

we have the two-point function of A p -A v at one-loop level as 

nf (po,tf T) = U^( Po ,p) + U i l- loop ^( Po ,p;T) . (3.7) 

Analogously to Eq. (|3.4|) . we split the two-point function into two parts as 

(p ,p;T) = U^(po,p)+U^( Po ,p;T) , (3.8) 

where n^ ac ^ l/ includes the quantum correction and the contribution at tree level in Eq. (|3.6|) . 
and represents the hadronic thermal correction. Since the hadronic thermal correction 
S^" has no divergences, the renormalization conditions for F 2 and zi can be determined from 
F 2 we adopt the "on-shell" renormalization condition: 

ni vac) ^(p = 0,p = 0) = g^F 2 (0) . (3.9) 

From this renormalization condition, we obtain the g^-part of rj^ ac ^ iy in the form |2j 

p, P ^r c ^( P0 ,p) = P 2 [f 2 (o) + ni( P 2 )] , (3.io) 

where 11^ (p 2 ) is the finite renormalization contribution satisfying 

Ul(p 2 = 0) = . (3.11) 
For Z2 we adopt the renormalization condition that TT^ ac - )At!/ be given by 



nf c) ^(p ,P) = <T [F'm +p 2 Tli(p I )\ + (g^p z -p» P n [2z 2 (M p ) +Ef (p*)j , (3.12) 

where Z2(M p ) is renormalized at the scale M p and Wf" {p 2 ) is the finite renormalization subject 
to the condition 

Refii T (p 2 = M 2 ) = . (3.13) 

To distinguish the hadronic thermal correction to the pion decay constant from that to the 
parameter we decompose the two-point function of A^-A v into four components as 

nf = uvni + - u ^u u )u s ± + pfTii + p^nl , (3.14) 
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Figure 2: Feynman diagrams contributing to the V M -V y two-point function. Here tt represents 
the quantum pion field and likewise for the others. 

where PP" and P^ 1 are the polarization tensors defined by 



{jgT - u»u a ) f -g af3 - t^L \ (gP» - U /V) 



^ ~ ) Ua ]pf u e [ 9 ~~^~J ■ (3 ' 15) 

Similarly to the division in Eq. (|3.8|) . it is convenient to divide each component into two parts 



as 



tf ± (Po,P;T) = nf ac)t ( Po ,p) + ni( P0! p;T) , (3.16) 

where Tl^ mc ^ t (po,p) includes the tree contribution plus the finite renormalization effect and 
n? ± (po,p;T) is the hadronic thermal contribution. With Eq. (|3,12|) the functions Yl^ ac " >t ' s ' L ' T 
can be written as 



r( vac ) L ^„ ^ _ rr( vac ) T /^„ x\ - _J2 



(P0,P) = Hi™ ; (Pd,p) = -P 2 [2. 2 (M p ) + 11^ (p*)j . (3.17) 

The explicit forms of the hadronic thermal corrections Tl t ^' L ' T (po,p; T) are summarized in 
Eqs. dEDJ-CUl) in AppendixEl 

We show the diagrams contributing to the V^-V u two-point function in Fig. [2j We adopt 
the on-shell renormalization condition similar to the above to obtain the resultant forms of four 
components of the two-point function as 
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n[ vac) Vp) = n[ vac) > ,p) = fI{m p ) + n£( P 2 ) 



n 



(vac) L , 



(po,p) = nS vac)T (p ,p) = -v l 2z 1 (m„) + nf V) 



(3.18) 



where the parameters are renormalized at the scale M p and the finite renormalization terms 
satisfy 



ReU^(p 2 = M 2 ) = ReUF(p 2 = M 2 ) = . 



(3.19) 



For the two-point functions of V P -V u and V^-Vu we adopt similar on-shell renormalization 
conditions. The resultant sums of the tree contributions and quantum corrections take the 
forms 



f 2 (m p ) + u s v (p 2 ) = nr i (Po,p) = n}^ fl (po,p) , 



r(vac)t 



(vac) s , 



ni; ac)L (po,p)=n(; ac)T (p ,p) 



+ n£ T (p* 



<7 )l (po,p) = (po,p) = V [«,(m,) + n-|(^) 

where, as in Eq. (|3.19|) . the finite renormalization terms satisfy 



r (vac)T 



(3.20) 



ReU^(p 2 = M l ) = ReU^ r (p 2 = M l ) = RsTL^W = M 



\LTf-2 



\LT/J2 







(3.21) 



The hadronic thermal corrections to the above two-point functions relevant to the present anal- 
ysis are given in Eqs. ()B.5|) and (|B.6|) in Appendix iBl 

It should be noticed that the renormalized parameters have the intrinsic temperature de- 
pendences in addition to the dependence on the renormalization point. Then, the notations 
used above for the parameters renormalized at on-shell should be understood as the following 
abbreviated notations: 



F v (0) = F n (n = 0;T) , 
F^Mp) = F^n = M p (T);T) , 
g(M p ) = g(p = M p (T);T) , 
z 1)2 ,3(M p ) = z 1)2j3 ^ = M p (T);T) , 



(3.22) 



where /i is the renormalization point and the mass parameter M p is determined from the on-shell 
condition: 
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Ml = M 2 JT) = <? 2 (/x = M p (T);T)F 2 (p = MJT); T) . (3.23) 



In addition, the parameter a appearing in several expressions in Appendices are denned as 

_ F2((i = M p (T);T) 



F^ = MJT);T) 



(3.24) 



4 Current Correlators 

We now turn to construct the axial-vector and vector current correlators from the two-point 
functions calculated in the previous section. The correlators are defined by 

G%(po = iu n ,p;T)6 ab = fJ T dr J d s xe~^ + ^ (j&(r, x) J£,(0, 0))^ , 

ip = iu n ,p; T)5 ab = dr J d^xe^^ ( J£(r, x)J%(0, 0))^ , (4.1) 

where J$ a and J£ are, respectively, the axial-vector and vector currents, iv n = 2nirT is the 
Matsubara frequency, (a, b) = l,...,Nj — 1 denotes the flavor index and ( )p the thermal 
average. The correlators for Minkowski momentum are obtained by the analytic continuation 
of p . 

For constructing the axial-vector current correlator (po,p;T) from the A p -A u two-point 
function, it is convenient to take the unitary gauge of the background HLS and parameterize 
the background fields £l and £r as 

= e"* , £r = e* , 4>= ~4> a T a ■ (4.2) 

where <j) denotes the background field corresponding to the pion field. In terms of this <p, the 
background A^ is expanded as 

A p = A p + d^ + --- , (4.3) 

where the ellipses stand for the terms that include two or more fields. Then, the axial-vector 
current correlator is 

GT = PaPp ^ +nf , (4.4) 

where the first term comes from the (^-exchange and the second term from the direct Ap-Ay 
interaction. By using the decomposition in Eq. (|3.14j) . this can be rewritten as 
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(4.5) 



where 



p z w ± u s ± 
-[p§n* ± -p2ni] 

-ifi + nT . 



+ n 



(4.6) 
(4.7) 



One can see from Eqs. (|4.6|) and (|4.7j) that the pion exchange contribution is included only in 
the longitudinal component G A . 

To obtain the vector current correlator Gy, we first consider the V propagator. By using 
the fact that the inverse V propagator i(D^ 1 )^ l/ is equal to U y v , the propagator for the field V 
can be expressed as 



- iD v v = u»u v D v + (ff 



'"' - u»u v )D v + PjfDp + P$ U D 



T 

V i 



(4.8) 



where 



P 2 (nf 



D 



V 



P 2 u v (u v - u v ) - p 2 u v (u v - Uy) ' 

P^y-Iiy) 



D 



v 



pllly{Wy - Uy) ~ p 2 Uy(Uy ~ Uy) ' 

w 



Dl=Dy 



p\Uy{Uy ~ Uy) ~ p^Uy (Uy ~ Uy) 
1 



U s - U T ' 
Ll v LL V 



(4.9) 
(4.10) 
(4.11) 
(4.12) 



By using the above propagator Dy and two-point functions of V^-V^ and V^-Vu, Gy can be 
put into the form 



Gy" =Uy^iDy ap Uy\+U^ . 



(4.13) 



After a lengthy calculation, we obtain 



G 



flV 

V 



n v 



p 

+ (g» v - u^u u ) 

W 

_ y H 

p2 



plUy^Uy - Uy) +p 2 (U v{] Uy ~ U V[] U 

i[l\^2 ll v\\ y- L v LL v\\ lL v L h 



-vl 



(UyUy^-Uy^Uy) + flL^ (lly ~ ify) ) } + U 
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^{-n£n y|| n y|| +n y|| (^v\\ u v + u h u v 



+ n 



- P 4K\\ (n^n-n^n) 



+ 



LL v\\ 
P 2 



(-pI (n y n y|| - n y|| n£) + p 2 n y|| (n y - n 



+ n 



T 



(4.14) 



One might worry that the above form does not satisfy the current conservation p^Gy U = 0. 
However since, as shown in Eq. ()B.5|) . the conditions 



TT* 



n 



V 



LL v\\ 
u s „ 



n{ 
nfi 



(4.15) 



are satisfied, Eq. (J4.14|) can be rewritten as 



D ^ (n y n y (n^ + 2n^ 



+ 



pIu v -p 2 u\ 



K„) ; 



n y (n£ + 211^) + (n 



+ n 



+ n 



(4.16) 



Now it is evident that the current is conserved since p^P 1 ^ = Pp,Pj" = 0. In the present analysis, 
the equality ILy = n y is seen to hold as shown in Eq. HB.5|) . This implies that II y = II y is also 
satisfied since the quantum corrections to II y and II y are equal to each other due to Lorentz 
invariance. Thus, Gy can be written as 



r L (jy -f- fj, (jy , 



(4.17) 



where 



G 



Gy 



Uy (Uy + 211 



LL v 



u L 



Uy (ill + 211 



+ n 



+ n 



Note that, in the above expressions, we have dropped the terms ^11 
are of higher order. 



and ^11 



(4.18) 
(4.19) 
since they 



14 



5 Pion Decay Constants and Pion Velocity 

We now proceed to study the on-shell structure of the pion. For this we look at the pole of 
the longitudinal component G\ in Eq. (|4.6|) . Since both 11^ and 11^ have imaginary parts, we 
choose to determine the pion energy E from the real part by solving the dispersion formula 







p2Reni(po,p;T)-p 2 Reni(p ,p;r)l (5.1) 



where p = \p\. As remarked in Section |21 in HLS at one-loop level, rPj_(po>P; T) and U s ± (po,p; T) 
are of the form 

niCp ,p;T) = ^(o) + ni( P 2 ) + ni( Po ,p;T) , 

ni(po,p;T) = F2(o) + ni(p 2 ) + ni(p ,p;T) , (5.2) 

where Hj_(p 2 ) is the finite renormalization contribution, and fPj_(j?o>£>; T) and IL^pojp^T) are 
the hadronic thermal contributions. Substituting Eq. 1)5.2(1 into Eq. (|5.1|) . we obtain 

0= (£ 2 -p 2 ) [F 2 (0) + Reni(p 2 = E 2 - p 2 )] + E^eU^E^; T) - p 2 Re tL s ± (E,p; T) . (5.3) 

The pion velocity v 7T (p) = E/p is then obtained by solving 



F 2 (0)+Reni(p,p;r) 
F2(0)+Refli(p,p;T) 



Here we replaced E by p in the hadronic thermal terms Y^^E^p) and II^_(£',p) as well as in 
the finite renormalization contribution 11^ (p 2 = E 2 —p 2 ), since the difference is of higher order. 
[Note that U'Kp 2 = 0) = 0.] 

Next we determine the wave function renormalization of the pion field, which relates the 
background field <j) to the pion field 7f in the momentum space as 

4> = Tt/F(p;T). (5.5) 

We follow the analysis in Ref. 23 to obtain 

F 2 (p;T) = Reni(£,;p;T) = F 2 (0) + Re fPj_ (p, p; T) . (5.6) 

Using this wave function renormalization and the velocity in Eq. (|5.4I) . we can rewrite the 
longitudinal part of the axial-vector current correlator as 

r L ( * p 2 U t ± (p ,p;T)U s ± (p ,p;T)/F 2 (p;T) L 

G A {p ,p) = ^ 2 2 — +U_ L {p ,p;T) , (5.7) 

- [Po ~ v 2 {p)p l + Tl n (p ,p; T)\ 
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where the pion self energy rL,-(po,p;T) is given by 
Il n (p ,p;T) 



1 (5-8) 



v\ {ni (p ,p; T) — Re ni (E, p; T) } - p 2 {III (p , P; T) - Re ni (£7, p; T) } 



Let us now define the pion decay constant. A natural procedure is to define the pion decay 
constant from the pole residue of the axial- vector current correlator. From Eq. (j5.7fl . the pion 
decay constant is given by 

Ht ± (E,p;T)W ± (E,p;T) 



LiP]T) F^p-T) 

[F 2 (0) + ni(p,p;r)] [F*(0) +U' ± (p,p;T)] ( _ () , 



F 2 (P;T) 

We now address how / 2 (p; T) is related to the temporal and spatial components of the pion decay 
constant introduced in Ref. |24j . Following their notation, let /* denote the decay constant 
associated with the temporal component of the axial-vector current and the one with the 
spatial component. In the present analysis, they can be read off from the coupling of the 7f field 
to the axial- vector external field A^: 

f t ( - T , _ tfAE,p;T) _ F 2 (0) + nj(p,p;r) 

fsfo T) = ^±(E,P;T) = F 2 (0) +nj (p,p;T) 
F{p-T) F(p;T) 

Comparing Eqs. (|5~TU|) and (|5~TT|) with Eqs. (fS"lj) . (f53|) and (JEHJ), we have [MI25| 

F(p;T) =Re/*(p;T) , (5.12) 
fl{p ] T) = fl{p;T)f°{p-T) , (5.13) 

* W Re/*(p;T) l ° J 

We are now ready to investigate what happens to the above quantities when the critical 
temperature T c is approached. Due to the VM in hot matter .4, the parametric p meson mass 
goes to zero (M p — > 0) and the parameter a approaches one (a — > 1), so we have [see Eq. ()B.7|) ] 

ni(p,p;T) T -. c -^j?(o ; r c ) = -§T 2 , 

fli(p,p;T) r -^— ^^(0;r c ) = -^T 2 . (5.15) 
Substituting these into the expression of the pion velocity in Eq. (|5.4|) . we obtain 
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w)^ 1 - ( 5 - 16 ) 

This is our first main result: in the framework of the VM, the pion velocity approaches 1 near 
the critical temperature * 6 , not as in the case of the pion-only situation |SJ[7j. 

From Eq. ()5.15|) . we can evaluate the pion decay constant Eq. (|5.9|) at the critical temperature 
which comes out to be 

AT, 

/2(p;T c ) = F2(0)-^T c 2 . (5.17) 

Since this /„- is the order parameter and should vanish at the critical temperature, the parameter 
F%(0) at T = T c is given at T c as 4 

AT. 

^» T ~^ Tc . (5.18) 

Substituting Eq. (|5.15j) together with Eq. (|5.18|) into Eqs. (|5.1U|) and (|5.11jl . we conclude that 
both temporal and spatial pion decay constants vanish at the critical temperature # 7 : 

/*(p;T c ) = /«(p;r c ) = 0. (5.19) 

This is our second main result. Note that while /*(T C ) = 0, xa(T c ) is non-zero in consistency 
with the lattice result. Here the HLS gauge boson plays an essential role. 



6 Axial- Vector and Vector Susceptibilities 



In terms of the quantities defined in the preceding sections, the axial-vector susceptibility 
Xa{T) and the vector susceptibility xv{T) for non-singlet currents * 8 are given by the 00- 
component of the axial-vector and vector current correlators in the static-low-momentum limit: 



Xa(T) 
Xv(T) 



2Nf lim lim 

1 p— 0p ->0 

2Nf lim lim 



G° v °(p ,p;T) 



(6.1) 



p— >0p — »0 

where we have included the normalization factor of 2Nf. Using the current correlators given in 



oo 



lim po ^ P 2 /P 2 



-1, we can express xa(T) 



Eqs. (|4.5j) and (|4.16j) and noting that lim po _ >0 P\ 

and xv (T) as 
* 6 Modulo small Lorentz-breaking correction mentioned above. 

* 7 These results differ from those obtained in a framework in which the only relevant degrees of freedom near 

chiral restoration are taken to be the pions [fJEl. We will explain how this comes about in the conclusion section. 
* 8 We will confine ourselves to non-singlet (that is, isovector) susceptibilities, so we won't specify the isospin 

structure from here on. 
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XA(T) 
Xv(T) 



2Nf lim lim 



p^Opo^O 



ni(po,p;T)-ni(po,p;r) 



2Nf lim lim 



p^Opo^O 



+ n 



(6.2) 



where for simplicity of notation, we have suppressed the argument (po,p;T) in the right-hand- 
side of the expression for \v{T). In HLS theory at one-loop level, the susceptibilities read 



X a(T) = 2N f 
XV (T) - 



2Nr lim lim 

1 p^Opo^O 



where the parameter o(0) is defined by 



F 2 (0) + lim lim {ni(po,p;r) - Ui( Po ,p;T)} 
p— +Opo^O l. J 

(a(0)F 2 (0) +w v ) (nfr + 2S 



$1 



a (o)F2(o) + n* v -n^ 



n 



(6.3) 



o(0) 



< vac)t (p = 0,p = 0) 



n 



(vac) f 
V 



( Po = 0,p = 0) 



(6.4) 



F%(0) F2(0) 
In Ref. |18j . a(0) was defined by the ratio F 2 (M p )/ F%(0) without taking into account the finite 
renormalization effect which depends on the details of the renormalization condition. In the 
present renormalization condition (|3.2(Jj) with Eq. (|3.21|) . the finite renormalization effect leads 
to 



U v {p 2 = 0) = -j^Ml (2 - v^tan" 1 y/s) 



and then a(0) reads 



o(0) 



F 2 AM P ) 



+ 



(4tt) 2 " 
N f Ml 



{2 - \/3tan _1 \/3 



(6.5) 



(6.6) 



F|(0) (470**2(0) 

It follows from the static-low- momentum limit of (flj_ — fi^) given in Eq. (|B.8|) that the axial- 
vector susceptibility xa(T) takes the form 



X a(T) = 2N f 



F 2 (0) - JV, J?(0; T) + iV> J 2 (M P ; T) 



-^/^{«/ 2 i(M p ;T)- P^OsT)} 



(6.7) 



Near the critical temperature (T — ► T c ), we have M p — > 0, a — ► 1 due to the intrinsic temperature 
dependence in the VM in hot matter [1]. Furthermore, from Eq. (|5.18j) . we see that the parameter 
F%(0) approaches ^fT 2 for T — > T c . Substituting these conditions into Eq. (|6.7|) and noting 



that 



lim 



X - 2 {J\{M P -T) - Jl^T)} 



M 2 
p 



-J?(0;T) 



24 ~ 



(6.8) 
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we obtain 

XA(T C ) = -fT 2 c . (6.9) 

To obtain the vector susceptibility near the critical temperature, we first consider a(0)-F 2 (0) + 
fly appearing in the numerator of the first term in the right-hand-side of Eq. ()6.3() . Using 
Eq. (jB.9|) . we get for the static-low-momentum limit of a(0)F%(0) + fly as 



lim lim 

p^Opo^O 



a(0)F^0)+%(po,p;T) 



Nf r ~ ~ ~ i 

= a(0)F 2 (0) - —j- [2J\(M P - T) - Jl{M p ; T) + a 2 J?(0; T)\ . (6.10) 

From Eq. (JEEJ) we can see that a(0) -> 1 as T ^ T c since F 2 (M p ) -> F 2 (0) and Af p 0. 
Furthermore, -F^(O) — ► -^T 2 as we have shown in Eq. (|5.18|) . Then, the first term of Eq. (|6.10|) 
approaches ^-T c 2 . The second term, on the other hand, approaches —-^T^ as M p — > and 
a — >• 1 for T — > T r . Thus, we have 



lim lim 

p^Opo— >0 



a(0)F 2 (0) +n* v (po,p; T)l — ♦ . (6.11) 



This implies that only the second term in the right-hand-side of Eq. 1)6.3(1 contributes to 
the vector susceptibility near the critical temperature. Thus, taking M p — > and a — * 1, in 
Eq. I)B.10(I . we obtain 

Xv(T c ) = -fT 2 c , (6.12) 

which agrees with the axial- vector susceptibility in Eq. 1)6.9(1 . This is a prediction, not an input 
condition, of the theory. For Nf = 2, we have 

Xa(T c ) = xv(T c ) = \Tl . (6.13) 

This is our third main result. The result Xv(T c ) = |T C 2 is consistent with the lattice result as 
interpreted in 9 . It is interesting to note that the RPA result obtained in [2j m NJL model 
in terms of a quasi-quark-quasi-antiquark bubble is reproduced quantitatively by the one-loop 
graphs in HLS with the VM. 

It should be noticed that the pion pole effect does not contribute to the ASUS in Eq. 1)6.91) 
since the pion decay constant /* vanishes at the critical temperature as we have shown in the 
previous section, and that the contribution to the ASUS comes from the non-pole contribution 
expressed in Fig. ^ In three diagrams, the third diagram in Fig. W^c) is proportional to (1 — a) 
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and then it vanishes at the critical temperature due to the VM. Similarly, since the transverse p 
decouples at the critical point in the VM PJE], the first diagram in Fig.^a) does not contribute. 
Thus, the above result for the ASUS in Eq. ()6.9|) comes from only the contribution generated 
via the (longitudinal) vector meson plus pion loop [see Fig. IHb)]. * 9 Similarly, the vector meson 
pole effect to the VSUS vanishes at the critical temperature as shown in Eq. ()6.11|) . and the 
contribution to the VSUS is generated via the pion loop [Fig. Hfc)] and the longitudinal vector 
meson loop [Fig. Hfb)]. * 10 Since the longitudinal vector meson becomes massless, degenerate 
with the pion as the chiral partner in the VM, loop contribution to the ASUS becomes identical 
to that to the VSUS. Thus, the massless vector meson predicted by the VM fixed point plays 
an essential role to obtain the above equality between the ASUS and the VSUS. 

In the present analysis, our aim is to show the qualitative structure of the ASUS and the 
VSUS in the VM, i.e., the equality between them is predicted by the VM. In order to compare 
our qualitative results with the lattice result, we need to go beyond the one-loop approximation. 
We note here that there is a result from a hard thermal loop calculation which gives xv(T c ) ~ 
1-3TJ? |22]. However this result cannot be compared to ours for two reasons. First we need to 
sum higher loops in our formalism which may be done in random phase approximation as in 
|26| . Second, the perturbative QCD with a hard thermal loop approximation may not be valid 
in the temperature regime we are considering. Even at T S> T c , the situation is not clear as 
pointed out in Ref. [2*7] . 

7 Summary and Remarks 

The notion of the vector manifestation in chiral symmetry requires that the zero-mass vector 
mesons be present at the chiral phase transition. As discussed by Brown and Rho the light 
vector mesons near the transition point "bottom-up" can be considered as Higgsed gluons in 
the sense of color-flavor locking in the broken chiral symmetry sector proposed by Berges and 
Wetterich |28l I29j and could figure in heavy-ion processes measured at RHIC energies. In this 
paper we are finding that in the VM, the vector mesons with vanishing masses at the chiral 

transition (in the chiral limit) can figure importantly in the vector and axial-vector susceptibili- 

* 9 Note that a in Fig. ^b) is the quantum field corresponding to the NG boson absorbed into the vector meson, 
i.e., the longitudinal vector meson. 

* 10 Note that the contribution from Fig. |5fa) vanishes since the transverse p decouples and that the one from 
Fig- Hid) also vanishes since it is proportional to (1 — a). 
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ties near the chiral transition point. The notable results are that the VM confirms explicitly the 
equality xv = XA at T c and that both /* and vanish simultaneously with the pion velocity 
v n ~ 1. These differ from the results expected in a scenario where only the pions are the relevant 
effective degrees of freedom. 

The reason for that the v n deviates from 1 is due to the Lorentz-breaking term in the bare 
HLS Lagrangian at which the matching to QCD is made in a thermal bath. We find the deviation 
is small (this will be detailed in [3U])- The small deviation from 1 is also found in dense skyrmion 
matter studied in 5 . 

If one assumes that the only light degrees of freedom near T c are the pions, then one can 
simply take the current algebra terms in the Lagrangian and the axial-vector susceptibility 
(ASUS) xa is uniquely given by the temporal component of the pion decay constant /* with the 
degrees of freedom that are integrated out renormalizing this constant. Then the unquestionable 
equality xv = XA at T c together with the lattice result Xv|t c leads to the Son-Stephanov 
result on the pion velocity v n = 0. There is however a caveat to this simple result and it is that 
the same reasoning fails when one computes explicitly the vector susceptibility (VSUS) using 
the same current algebra Lagrangian. 

Positing that the vector mesons enter in the VM near T c circumvents this caveat and at the 
same time, makes a concrete prediction. In this framework, the ASUS is given by a term related 
to f* plus contributions from the vector-meson (i.e., the longitudinal component a) loop. At 
T c , the f* vanishes and what remains comes out precisely equal to the VSUS xv i n which the a 
(longitudinal p) loop in xa is replaced by the pion loop. All these are perfectly understandable 
in terms of the VM in HLS. 

If chiral symmetry restoration a la HLS/VM - but not the standard chiral theory one - is 
the valid scenario as is argued in [2] , its confirmation would provide a valuable insight into some 
of the basic tenets of effective field theories as expounded in [2J. 

As was shown in, e.g., Ref. |1(J1 we expect that baryons become light very near the 
phase transition point. At least when matter density is involved, the light baryons must figure 
importantly. In |31j . we have simulated density effects by introducing quasiquarks whose effective 
mass is expected to lie below the in- medium vector meson mass m*. In [S], the in-medium vector 
mesons whose mass must be much higher than that of quasiquarks were integrated out and the 
SUS was then given by the RPA bubble of quasiquark-quasiantiquark excitations in NJL model. 
In this paper where the vector meson plays the key role, it is not clear that we are not double- 
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counting if we introduce both the vector meson and the quasiquark in HLS/VM theory. Our 
point of view here is that whatever fermionic degrees of freedom are to be implemented in 
the theory should be color-singlet and hence the relevant fermionic degrees of freedom must 
be baryons. Assuming that the baryons are always more massive in the heat bath, our results 
should then correspond to HLS/VM with the baryons integrated out. Though our results clearly 
indicate the dual nature of the fermionic RPA bubble |U] and the one-loop HLS/VM, how to 
consistently implement fermions in a calculation of the type we considered here is still an open 
issue. We intend to return to this in a later publication. 

So far we have not addressed the properties of the quantities we have studied in this paper 
away from the critical point T c . The "intrinsic dependence" crucial in our formalism is a difficult 
problem to solve away from the T = and T = T c points. We have not yet formulated how to 
go about this problem. Confrontation with future lattice data as well as with RHIC data will 
obviously require these properties to be worked out. In the case of density, this problem has 
been approached from a different perspective in [321 • ^ is plausible that a similar approach can 
be developed for the temperature case. 

A more comprehensive discussion of the materials covered in this paper as well as other issues 
of HLS-VM in hot bath near chiral restoration will be discussed in a future publication |14j . 

In the present analysis, we used the bare HLS Largangian with Lorentz invariance since 
the Wilsonian matching between the HLS and the underlying QCD showed that the Lorentz 
violating effects to the bare parameters of the HLS Lagrangian are small, as we discussed around 
the end of section [2] and Appendix EI The details of the inclusion of such small corrections at 
the bare level and quantum effects based on such a Lagrangian will be presented in future 
publications. 
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Appendices 



A Wilsonian Matching and Intrisic Temperature Dependences 

The Wilsonian matching was originally proposed at T = 18 to determine the bare parameters 
of the HLS by matching the HLS to the underlying QCD. In Ref. [I] the Wilsonian matching was 
extended to non-zero temperature and it was shown that the parameters of the HLS Lagrangian 
have the intrinsic temperature dependences. In this appendix we first briefly review the Wilsonian 
matching proposed in Ref. ^H] at T = to determine the bare parameters of the HLS Lagrangian 
by matching the HLS with the underlying QCD. (For details, see Ref. 0.) Then, we extend 
the Wilsonian matching to the analysis at non-zero temperature to determine the intrinsic 
temperature dependences of the bare parameters of the (bare) HLS Lagrangian needed in the 
present paper. 

A.l Wilsonian matching conditions at T = 



The Wilsonian matching proposed in Ref. ^E] is done by matching the axial-vector and vector 
current correlators derived from the HLS with those by the operator product expansion (OPE) 
in QCD at the matching scale A. * 11 The axial-vector and vector current correlators in the 
OPE up until 0(1/Q e ) at T = are expressed as [33] 



nr- m = ^ (f ) [- ( 



8iV c vr ) /x 2 



| ^ {^G, u Gn | ^ 96(iV 2 - 1) n | l\a s (qq) 



N c Q 4 N c N% \2 3N C J Q 



(A.l) 



r (QCD),^ 1 fN r 



% W ) - 8?r 2 I 3 



3(^-l)M ln ^ 
SN C vr J ^ 



+ vr 2 (^G^Gn vr 3 96(iV 2 - 1) (I a s (qq) 



(A.2) 



N c Q 4 N c N% \2 3N C J Q e 

where /i is the renormalization scale of QCD and we wrote the ^-dependences explicitly (see, 
e.g., Ref. 34 ). In the HLS the same correlators are well described by the tree contributions 
with including 0(p 4 ) terms when the momentum is around the matching scale, Q 2 ~ A 2 : 

n? LS) (Q 2 ) = ^!^-2z 2 (A) , (A.3) 



* n For the validity of the expansion in the HLS the matching scale A must be smaller than the chiral symmetry 
breaking scale A x . 
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n 



(HLS) 



(Q 2 



^ 2 (A) 

m 2 (A) + q 2 



l-2g 2 (A)z 3 (A) -2«i(A) 



(A.4) 



where we defined the bare p mass M p (A) as 



M 2 (A) ^ 5 2 (A)F 2 (A) . 



(A.5) 



We require that current correlators in the HLS in Eqs. (|A.3|) and (|A.4|) can be matched with 
those in QCD in Eqs. (fA~T|l and (TOl . Of course, this matching cannot be made for any value of 
Q 2 , since the Q 2 -dependences of the current correlators in the HLS are completely different from 
those in the OPE: In the HLS the derivative expansion (in positive power of Q) is used, and the 
expressions for the current correlators are valid in the low energy region. The OPE, on the other 
hand, is an asymptotic expansion (in negative power of Q), and it is valid in the high energy 
region. Since we calculate the current correlators in the HLS including the first non-leading 
order [0(p 4 )], we expect that we can match the correlators with those in the OPE up until the 
first derivative. Then we obtain the following Wilsonian matching conditions |181 [2] * 12 



^ 2 (A) 
A 2 



1 

8^ 



N, 



1 + 



3(iV 2 -l) a £ + 2vr 2 <^GH 



+ 



8N C 
288tt(N 2 - 1 



7T 



N2 



N c 
1 



1 

2 + 3N C 



A 4 

u s (qqY 
A 6 



(A.6) 



F 2 (A) A 4 [l-2g 2 (A)z 3 (A)] 
(M p 2 (A) + A 2 ) 2 



A 2 



_L (El 

8tt 2 { 3 



3(JV C 2 - 1) a s 2tt 2 (^G^CT) 

8N r 7T N r A 4 



288vr(A r , 



1) 



AT3 

F 2 (A) F 2 (A)[l-2g 2 (A)z 3 (A)] 
A 2 M p 2 (A) + A 2 

4tt(N 2 - 1) a s (qq) 2 



1 



3N C 



u s (qgY 
A 6 



(A.7) 



2[z 2 (A)-Z!(A)] 



(Ai 



A 6 " 

The above three equations l|A.8|l . (|A.6|) and ()A.7|) are the Wilsonian matching conditions pro- 
posed in Ref. ^S] . These determine several bare parameters of the HLS without much ambiguity. 

Especially, the first condition HA.6j) determines the ratio i ? vr (A)/A directly from QCD. 
* 12 One might think that there appear corrections from p and/or 7r loops in the left-hand-sides of Eqs. 1A.6I and 
(S3- However, such corrections are of higher order in the present counting scheme, and thus we neglect them 
here at Q 2 ~ A 2 . In the low-energy scale we incorporate the loop effects into the correlators. 
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A. 2 "Intrinsic" temperature dependence of the bare parameters 

Let us consider the extension of the above matching conditions to the analysis in hot matter. The 
quark condensate as well as the gluon condensate appearing in the right-hand-side (RHS) of the 
above matching conditions generally have the temperature dependence which is converted into 
the intrinsic temperature dependence of bare parameters through the matching conditions [I]. 
We should note that there is no longer Lorentz symmetry in hot matter, and the Lorentz non- 
scalar operators such as q^y^Dyq may exist in the form of the current correlators derived by the 
OPE ^]. Such a contribution may generate Lorentz symmetry violating effects in the RHS of 
the above matching conditions, and accordingly, we may have to use the bare Lagrangian with 
the Lorentz symmetry breaking effects included as in Appendix A of Ref. |31| . However, we 
will see that we can use the Lorentz invariant form of the bare Lagrangian in Eq. ()2.5j) near the 
critical temperature as a good approximation as follows: 

In the RHS of the matching condition in Eq. (|A.6|) . the Lorentz symmetry violating contri- 
bution from the operators such as q , y l j i D u q are small compared with the main term of 1 + 
This implies that the Lorentz symmetry breaking effect in the left-hand-side of Eq. (|A.6|) , which 
is expressed by the Lorentz symmetry violation in the bare ir decay constant, is also small: The 
difference between F* bare and F£ bare is small compared with their own values, or equivalently, 
the bare vr velocity defined by v\ bare = F° hai jF^ 

bare * s c l° se to one. As a result we can deter- 
mine, in a good approximation, the bare ir decay constant through the matching condition in 
Eq. I)A.6|) with putting possible temperature dependence on the gluon and quark condensates [I]: 



where we took N c = 3. We should stress again that, through the above condition (|A.9|) . the tem- 
perature dependence of the quark and gluon condensates determines the intrinsic temperature 
dependence of the bare ir decay constant F 7T (A;T). 

Next, we consider the intrinsic temperature dependence of other parameters near the critical 
temperature. As was shown in Ref. [HJ f° r the VM in dense matter, the equality between the 
vector and axial- vector current correlators in the HLS requires the following VM conditions for 
the bare parameters at the leading order, which should be valid also in hot matter: 



^ 2 (A;T) 
A 2 




a s 2vr 2 {^G^Gnr 

7T 3 A 4 



+ 7T 



3 1408 a s {qq) 2 T 
27 A 6 



(A.9) 




(A.10) 



5T,bare 







9 L, bare 







for T -^T c . 



(A.ll) 
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The VM conditions for the a parameter in Eq. (jA. 10|) together with the above result that the 
Lorentz symmetry violating effect between F* bare and bare is small already implies that the 
effect of Lorentz symmetry breaking between F* bare and F* bare is small: The bare velocity of a 
(longitudinal p) defined by v 2 haxe = F* bare /F* bare is close to one near the critical temperature 
determined from the intrinsic temperature dependence. On the other hand, the ratio t>T,bare = 
9 L, bare / 9t, bare, which we call the bare velocity of the transverse p, cannot be determined through 
the Wilsonian matching, since the transverse p decouples near the critical point in the VM ^ 
121 Hj. However, this decoupling nature of the transverse p near the critical temperature implies 
that it becomes irrelevant to the quantities studied in this paper. Thus, in the present analysis, 
we set ^T.bare = 1 for simplicity of the calculation, and show how the transverse p decouples 
from the quantities we study in this paper. 



B Hadronic Thermal Corrections 

In this appendix we summarize the hadronic thermal corrections to the two-point functions 
of A^-A u , V^-Vv, V^-Vu and F M -V„. 

The four components of the hadronic thermal corrections to the two point function of A p -A u , 
are expressed as 

ni(p ,p;T) = N f (a - 1)A (0,T) - N f aM 2 p B Q (p Q ,p- M p ,0;T) 

+ N f ^B t (p ,p;M p ,0;T) , (B.l) 

ni(po,p;T) = N f (a - l)I (0,r) - N f aM 2 p B (p ,p;M p ,0;T) 

+ N f ^B s (p ,p;M p ,0;T) , (B.2) 

ni(p ,p;T) = N f jB L (p ,P;M p ,0;T) , (B.3) 

tll(p ,p;T) = N f ~B T (p o ,p;M p ,0;T) , (B.4) 

where the functions Ao, Bq, and so on are given in Appendix IU1 

The two components IF and IT of hadronic thermal corrections to the two-point functions 
of V p -V„, Vfj,-V u and V p -V u are written as 

U t v (p ,p;T)=n s v (po,p;T) 
= n\( Po ,p;T) = Ul(po,P;T) 
= -& v ,,(po,p;T) = -tl s vll (p ,p;T) 
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1 



-N f - A (M p ;T)+a'A (0;T) - N f M z B (p ,p; M p , M p ;T) . 



(B.5) 



Among the remaining components only is relevant to the present analysis. This is given 
by # 13 



1 



If (p ,p; T) = N f -B L (p ,p; M p ,M p ;T) + N f 



(2-q) 

8 



-B L (p ,p;0,0;T) 



(B.6) 



For obtaining the pion decay constants and velocity in Section |S] we need the limit of po = p 
of fF ± and U s ± in Eqs. (|B.lj) and l|B.2j) . With Eq. (|C.6|) . and U s ± reduce to the following 



forms in the limit M p — ► and a —> 1: 



ni(p =p + ie,p;T) 
nl(p =p + ie,p;T) 



M p ->0, o->l 



-^J 2 (0;T) = -^T 2 , 
2 iv > ; 24 

2 iv i ; 24 



(B.7) 



In the static-low- momentum limits of the functions listed in Eq. (|C.7|) . the (IIj_ — 11^) 
appearing in the axial-vector susceptibility becomes 



lim lim 

p^Opo— >0 



ni(po,p;T)-ni(po,p;T) 



-N f jf(0;T) +N f a jf(M p ;T) - N fJ ^ [J^(M p ;T) - ^(OjT) 



(Bi 



For the functions appearing in the vector susceptibility relevant to the present analysis we have 



lim lim 

p— >Opo— >o 

lim lim 

p^Opo— >0 



n' y (p ,p;T) 
'njf(po,p;r) 



-Nt 



2J° 1 (M p ; T) - jf(M p ; T) + a 2 J 2 (0; T) 

2 



M 2 p 4(M p -T) + 2Jf{M p -T) 



(2 -a) 



-Ji(0;T) . 



(B.c 



(B.10) 



C Functions 

In this appendix we list the explicit forms of the functions that figure in the hadronic thermal 
corrections, Aq, Bq and B^ v in various limits relevant to the present analysis. 
The function Aq(M;T) is expressed as 

A (M;T) = J?(M;T) , (C.l) 

where J 2 (M;T) is defined by 

* 13 The explicit forms of other components will be listed in Ref. |14). 
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JP(M;T) 



d 3 k 



\k\ 



n-2 



(2vr) 3 e ^(fc;M)/T _ ! [ u (k;M)] 1 



(C.2) 



with I and n being integers and uj{k;M) = y M 2 + \k\ 2 . In the massless limit M = 0, the 
above integration can be performed analytically. Here we list the result relevant to the present 
analysis: 



(C.3) 



It is convenient to decompose B^ v into four components as done for 11^ in Eq. (|3.14|) : 



= u^B 1 + (g^ - u^u u )B s + Pj* u B + PfB 



(C.4) 



We note here that, by explicit computations, the following relations are satisfied: 



S*(po,p;M,M;T) = B s ( Po ,p;M,M;T) = -2A Q (M;T) = -2jf(M;T) . 



(C.5) 



To obtain the pion decay constants and velocity in Section |S] we need the limit of po = p 
of the functions in Eqs. ()B.1|) and (|B.2|) , As for the functions M 2 Bq, B l and B s appearing in 
Eqs. (|B.1|) and (|B.2|) . we find that, in the limit of M p going to zero, they reduce to 



M 2 B (p = P + ie, p; M p , 0; T) -— , 

B t (p =p + ie,p;M p ,0;T) — ► -2J?(0; T) = -\t 2 , 

B s (p =p + ie,p;M p ,0;T) — ► -2J?(0; T) = -V . 



(C.6) 



The static-low-momentum limits of the functions appearing in the corrections to the axial- 
vector and vector susceptibility are summarized as 



lim lim 

p— >0po^0 

lim lim 



M 2 p B (p ,p; M p , 0; T) = -j\{M p - T) + J 2 (0; T) 



B'ipcP; M p , 0; T) - B L (p ,p; M p , 0; T) 



-4 r 



M 2 



-J 2 _ 1 {M p ;T) + Jl^-T) 



lim lim 

p— >Opo— >o 

lim lim 

lim lim 

p^Opo^O 



M Z pB (p ,p;M p ,Mp;T) =- J^(M p ;T) - J{(M p ;T) 



B L ( Po ,p;Mp,M p ;T)\ = -2M 2 p J?(M p ;T) - 4Jf(M p ;T) , 
B L (p ,p;0,0;T)} = -AJ 2 (0;T) . 



(C.7) 
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